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Fig. 1. Magnetic fields and spin periods of observed pulsars (data 
taken from the ATNF pulsar catalogue). Black dots are normal 
pulsars. Red dots represent MSPs. Green triangles are pulsars 
in binaries. The "spin-up line" represents the minimum spin- 
period to which a spin-up process may proceed in an Eddington- 
limited accretion, while the "death-line" corresponds to a polar 
cap voltage below which the pulsar activity is likely to switch 
off (Bhattacharya & van den Heuvel 1991). 



Radio pulsars are found to be of two kinds, i.e., normal pul- 
sars, with magnetic field B ~ lO'^ G and spin period P ~ a few 
seconds, and millisecond pulsars (MSPs), which have low mag- 
netic fields {B ~ 10^ - 10'^ G) and short spin periods (P < 20 
ms, e.g. Bhattacharya & van den Heuvel 1991; Lorimer 2008). 
These bimodal distributions in magnetic fields and spin periods 
are shown, e.g., in Fig. 1 in Wang et al. (2011) and also in FigH] 
of this paper, which, with data taken from ATNF pulsar cata- 
logue, shows a large population of normal pulsars and a smaller 
population of MSPs. These two populations are connected with 
a thin bridge of pulsars in binaries. Moreover, most MSPs are 
in binary systems. It has long been proposed that MSPs are 
formed through a recycling process in which neutron stars ac- 
crete material from their low-mass companions and are spun 
up by the angular momentum carried by the accreted material. 
During this phase, the magnetic field is buried by the accreted 
material and decays (Alpar et al. 1982; Taam & van den Heuvel 
1986; Bhattacharya & van den Heuvel 1991; Radhakrishnan 
& Srinivasan 1982; Bhattacharya & Srinivasan 1995; van den 
Heuvel 2004). These objects are therefore called recycled MSPs 
(Taam & van den Heuvel 1986; Bhattacharya & van den Heuvel 
1991; van den Heuvel 2004). 

The evolutionary precursors of recycled MSPs are widely 
believed to be pulsars in binaries with high magnetic field 
(B ~ 10" - 10'^ G) and long spin period (P ~ a few sec- 
onds, Bhattacharya & van den Heuvel 1991; Liu, van Paradijs 
& van den Heuvel 2007; Lorimer 2010). It is evident that B 
and P of normal pulsars and recycled pulsars are correlated with 
the duration of accretion phase and the total amount of accreted 
mass (Taam & van den Heuvel 1986; Shibazaki et al. 1989; 
Wijers 1997). If a neutron star accretes a small amount of mass 
from its companion, e.g. ~ O.OOIMq - O.OIMq, a recycled pul- 
sar with mildly weak field and short spin period {B ~ 10'" G, 
P ~ 50 ms) will be formed (e.g. Francischelli, Wijers & Brown 
2002), like PSR 1913+16 and PSR J0737-3039 (Lyne et al. 
2004; Kramer 2006). Strong supporting evidences for this recy- 
cling idea have been found in low mass X-ray binaries (LMXBs) 



containing accreting millisecond X-ray pulsars (AMXPs), e.g. 
SAX J 1808.4-3658 (Wijnands & van der KHs 1998), and in 
observing the transition link from an X-ray binary to a radio 
pulsar PSR J1023+0038 (Archibald et al. 2009). At the end of 
the accretion phase (accreted mass greater than 0.2Mq), the neu- 
tron star magnetic field may arrive at a bottom value of about 
10'^ - 10^ G and its spin period may reach a minimum of about 
a few milliseconds. A millisecond pulsar is formed (van den 
Heuvel & Bitzaraki 1995ab; Urpin, Geppert & Konenkov 1998). 
The accretion-induced field-decay and spin-up torque make the 
magnetic field and spin period (Hobbs & Manchester 2004; 
Manchester et al. 2005) to change from B ~ 10" - 10"'' G and 
P ~ a few seconds to B ~ 10** - 10'' G and P ~ a few mil- 
liseconds. It is the recycling process that leads to the bimodal 
distribution of radio pulsars. 

In this paper, we examine the accretion-induced field-decay 
and spin-up model of Zhang & Kojima (2006) for these recy- 
cled pulsars. We investigate the differences in model-predicted 
B and P distributions with different initial conditions. We also 
compare those distributions with currently observed ones. This 
paper is organized as the following. Section 2 gives an overview 
of the model. We describe the input parameters necessary in the 
calculations of field and spin evolutions and analyze the results 
in Section 3. Section 4 contains discussions and summary. 

1. The model 

1.1. The recycling process 

A neutron star in a binary, with an initial magnetic field of about 
Bo ~ 10'^ G and an initial spin period of about Pq ~ a few sec- 
onds, can accrete material from its companion and forms an ac- 
cretion disk. When the ram pressure of accretion material equals 
the magnetic pressure, a magnetosphere forms. As a result, a 
boundary layer appears between the innermost disk and mag- 
netosphere due to the transition for rotating velocity of plasma 
from Keplerian to the spin velocity of the neutron star (Inogamov 
& Sunyaev 1999). In this layer, the accreted matter will be chan- 
neled onto the polar patches by the field lines, where the com- 
pressed accreted matter causes the expansion of magnetic polar 
zone in two directions, downward and equatorward (Zhang & 
Kojima 2006). Therefore, the magnetic flux in the polar zone is 
diluted, and more matter is accreted to the polar cap and diffuses 
to the surface of the neutron star. Finally, the polar cap area ex- 
pands and occupies the entire neutron stars surface, and the mag- 
netic flux is buried in the equatorial area. The magnetosphere is 
then compressed to the neutron star surface, leading to an object 
with weak fields in large scale (about ~ 10** G) and very strong 
fields in small scale (about ~ 10'^ G). Meanwhile, the angular 
momentum carried by the accreted matter spins up the neutron 
star, forming a MSP with the spin period of a few milliseconds. 



1.2. Magnetic field evolution 

Based on the above accretion-induced field-decay and spin-up 
model, the accretion-induced field and spin evolution is obtained 
analytically: (Zhang & Kojima 2006) 
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xl - (l^)"*^^- Bq - B(t - 0) is the initial field sti-ength and 
Bf is the bottom magnetic field, which is defined by the neu- 
tron star magnetosphere radius matching the stellar radius, i.e., 
RM(Bf) = R. Rm is defined as Rm = (pRA where Alfven radius 

Ra = 3.2 X lO'^cmMj^^^Vto^'w"'^^ (Eisner & Lamb 1977; Ghosh 
& Lamb 1977). The model dependent parameter ^ is about 0.5 
(Ghosh & Lamb 1979b; Shapiro & Teukolsky 1983; Frank et al. 
2002). Mi7 is the accretion rate in units of 10'^ g/s. //so is the 
magnetic moment in units of lO-'" G cm-*. The mass m - M/Mq 
is in the unit of solar mass. According to this model, the bot- 
tom field of neutron stars is determined by the condition that the 
magnetosphere radius equals the neutron stars radius (Zhang & 
Kojima 2006). Using the relation RM{Bf) = R, we can obtain the 
bottom field, 

Bf ^ 1.32 xlQ^G(^)^- m-*R'J(f>-\ (2) 
Mi 8 

where Mi8 = M/IO"* g/s und R(, = R/IO^ cm. 



1.3. Spin evolution 

During the accretion phase, the neutron star is spun up by the 
angular momentum carried by the accreted matter. The spin 
evolves according to the following relation given by Gosh & 
Lamb (1979b): 

-P = 5.8xlO-'[(^r^Rjl,^] 



(3) 



where we define the parameters, the surface field B\2 - B/IO'^ 
G, the moment of inertia 745 - //lO'^^g cm^, the X-ray bright- 
ness (L - GMM/R) L37 in units of lO"*^ erg/s, respectively. The 
dimensionless parameter n(wj) is the fastness parameter, whose 
expression is given by Gosh & Lamb (1979b), 
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with Q.S being the stellar spin frequency and Q.k is the Keplarian 
frequency. Wj is the ratio parameter of the angular velocities 
which describes the relative importance of stellar rotation and 
plays a significant role in our entire understanding of accretion to 
the rotating magnetic neutron stars (Eisner & Lamb 1977; Ghosh 
& Lamb 1977; Li & Wang 1996, 1999; Shapiro & Teukolsky 
1983). For a slowly rotating magnetic neutron star, <k 1. 
Uc depends on several properties of the system (Gosh & Lamb 
1979b) and is taken to be 0.35 in this computation. 

2. B and P distributions of IVISPs 

2.1. Input parameters 

According to the accretion-induced field-decay and spin-up 
model by Zhang & Kojima (2006), a pulsar with strong magnetic 
field (e.g. B ~ 10'^ G) and slow rotation {P ~ a few seconds) in 
a binary system may be spun up to become a millisecond pulsar 
via accretion, and the magnetic field decays to a bottom value 
(B ~ 10** - 10'' G) during this phase. The final state of a recycled 



system is characterized by the magnetic field and spin period, 
which are related to the initial magnetic field, initial spin period, 
accretion rate, accretion time, and the mass and radius of neutron 
stars. In computing the model-prediction of B and P, we adopt 
these input parameters as the following: 

(1) All precursors to recycled MSPs are assumed to be normal 
pulsars with magnetic field of about Bo ~ 10'^ G and spin period 
of about a few seconds. According to the data taken from ATNF 
catalogue, the log B and log P distributions of normal pulsars 
can be well described by a Gaussian function (see Wang et al. 
201 1). We therefore take lognormal distributions as the inputs of 
initial magnetic fields and spin periods. The probability density 
function reads. 



p(x) 
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where fi is the mean value of the distribution, and cr is the stan- 
dard deviation. We take fj. = logB(G) = 12 and cr - 0.5 for 
the B distribution (see e.g. Hartman et al. 1997; Hobbs et al. 
2011; Kaspi 2010; Wang et al. 2011) and = logP(s) = and 
cr = 0.4 for the P distribution (see e.g. Lorimer 2010 and ref- 
erences therein). The range of B and P to consider is taken as 
Bo = lO'"-^ - 10'4 " G and Pq = 0.1 - 30 s, respectively. 

(2) The accretion rate is taken in the range from M - 10'^ g/s 
to M = 10'** g/s (see e.g. Wijers 1997) and is assumed to be 
constant during the whole process. The duration of the accretion 
phase is in the range from A? = 10^ yr to Af = 10'^ yr. Most sys- 
tems may accrete 0.1- 0.2Mo at the end of the accretion phase 
(see Shapiro & Teukolsky 1983). However, in some binary sys- 
tems, the maximum accretion mass can be 0.6 - O.8M0 (van den 
Heuvel 2011), depending on the intrinsic properties of the sys- 
tem. We adopt lognormal distributions (see Eq.|6]l for the accre- 
tion rate and the accretion time. The mean values and standard 
deviations are fi = logM(g/s) - 17, cr = 0.1 for the accretion 
rate and fi - log Af(yr) = 8, cr = 0.4 for the accretion time re- 
spectively. 

(3) According to recent statistics of neutron star mass (see 
Zhang et al. 2011), we consider a Gaussian mass distribution 
with the mean at 1 .4M0 and the standard deviation equal to 
0.2 Mq within the range from 0.9 Mq to 2.2 Mq. 

(4) It is widely believed that the neutron star radius is about 10 
km. We consider a uniform distribution of the radius from 10 km 
to 20 km. 

2.2. Distribution of recycled MSPs 

We calculated the final B and P according to the accretion- 
induced B and P evolution model described in section [T] In or- 
der to produce the model-predicted distribution, we generated a 
large number of input parameter random samples from the in- 
put parameter distributions mentioned in section 12.11 Different 
computation runs with 300, 3000 and 30000 samples were per- 
formed to examine how many samples are needed to have a con- 
verged result. Fig.|2]shows the B and P distributions of recycled 
MSPs obtained with 300, 3000 and 30000 random samples, re- 
spectively. We can see that, with the initial conditions consid- 
ered in section ITTl the Zhang & Kojima (2006) model indeed 
expects the magnetic fields of recycled MSPs are of 10^ - 10'^ G, 
and their spin periods are a few milliseconds. We divided these 
distributions into histograms with bin size of 0. 1 and compared 
their difference with the Kolmogorov-Smirnov (K-S) test (Press 
et al. 1992). For the comparison of B distributions between 300 
and 3000 samples, the K-S statistic gives a 86.1% probability 
of the two being drawn from the same parent distribution. That 
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Fig. 3. B and P input distributions with narrow and wide Gaussian widths (the upper panels) and their resuhant final B and P 
distributions (the lower panels). 
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Fig. 2. The computed B and P distributions of recycled MSPs 
with different numbers of random input samples. The left (right) 
panel is the B iP) distribution. As discussed in the text, the com- 
puted distributions with 30000 random samples are good enough 
and are the ones used for comparison with observation in this 
work. 



usually defined as pulsars with spin period less than 20 ms. Fig. 
|4] shows the B and P distributions of observed MSPs. Because 
of the wide ranges of the input parameters, we have some MSPs 
with longer computed spin periods (see Fig.|2|. The comparison 
with observation is conducted for spin periods shorter than 20 
ms, and the computed B and P distributions for these recycled 
pulsars are plotted in Fig.|5] Again using the K-S test, we found 
that the degree of similarity for the B distributions is 69.2% and 
that of the P distributions is 73.6%. We also performed ?lx^ test 
to check the consistency between the model predicted distribu- 
tions and observed ones. The errors in the observed distributions 
are assigned as V^V; for each bin, where A^, is the number of ob- 
served recycled pulsars in the /th bin. The;^f^ value is 22.8 (19 
degrees of freedom) for the B distribution, which corresponds 
to a null hypothesis probability of 24.6%, and 16.6 (11 degrees 
of freedom) for the P distribution, which corresponds to a 12.0% 
null hypothesis probability. Both the K-S test and;^^^ test indicate 
that the consistency between the currently observed distributions 
and the Zhang & Kojima (2006) model prediction is roughly at 
about the so-called 1-cr level. In other words, although not of a 
high degree of similarity, they are not inconsistent to each other 



probability for the two distributions of 3000 and 30000 samples 
is 100.0%. For the P distributions, that probabiHty is 64.9% for 
that of 300 and 3000 samples, and again reaches 100.0% for that 
of 3000 and 30000 samples. We therefore content ourselves with 
the results obtained with 30000 random samples. 

Although the initial conditions, that is, input parameters, as 
described in section 12.11 are quite commonly adopted, we fur- 
ther checked the sensitivity of the model prediction results to 
the variation in the distribution of the input parameters. We used 
different input widths of the lognormal B and P distributions, 
which are 0.375 (instead of 0.5) for B and 0.2 (instead of 0.4) 
for P. The ranges exploited for the B and P distributions are the 
same as in section ITTI These input distributions and the corre- 
sponding computed final B and P distributions are shown Fig. [3] 
Based on the K-S test, the similarity between the final distribu- 
tions obtained with narrow and wide inputs is 100.0%, for both 
B and P distributions. They are insensitive to mild variation of 
their progenitor distributions. 

Finally, we compared the model-predicted B and P distri- 
butions with currently observed ones. The observed MSPs are 



2.3. B and P evolutions of recycled MSPs 

Fig. |6] displays the B and P evolutionary tracks during the ac- 
cretion phase with initially slowly rotating and highly magne- 
tized progenitors. In this calculation, we took the initial mag- 
netic field as Bo = 5 X 10'^ G and the initial spin period as 
f() = 1 s. Accretion rates of 10"* g/s, 10'^ g/s, and 10'^ g/s were 
considered. These tracks were followed until the accreted mass 
reached 1 Mq. The upper right panel in Fig. |6] shows that the 
magnetic field decays with the accumulation of the accreted ma- 
terial and that the bottom values of magnetic fields correlate with 
the accretion rates. The lower left panel shows that the pulsar ro- 
tates faster and faster when accreting more and more mass and 
the spin period is insensitive to the accretion rate after accreting 
about O.OOIMq. 

3. Discussion and Summary 

We tested the accretion-induced field-decay and spin-up model 
for recycled MSPs (Zhang & Kojima 2006, Wang et al. 2011). 
In our computation, we considered lognormal distributions of 
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Fig. 6. B and P evolution in the recycling process. The upper left panel shows the joint evolution of B and P. The upper right and 

lower left panels are their evolution as a function of accreted mass AM. The solid, dashed, and dotted lines are the evolutionary 
tracks with the accretion rate of 10"^ g/s, 10^^ g/s, and 10'* g/s, respectively. The initial B and P are taken as Bo = 5 x 10'^ G and 
Po = 1 s. The lower right panel is a zoom-in view of the lower left panel in a Unear scale for spin periods shorter than 5 ms. 



initial and P in the range of Bq = lO'^-^ - lO'^ G and = 
0.1 - 30 s, a Gaussian neutron star mass distribution in the range 
of M = O.9-2.2M0, a uniform neutron star radius distribution in 
R = 1 - 20 km, , a lognormal distribution of the accretion rate in 
M = 10'* - 10'^ g/s and of the accretion time in At = 10^ - lO' 
yr. We found that the computed B and P distributions of recycled 
MSPs are insensitive to mild variations in the width of the initial 
distributions. Based on the K-S test and test, we found that 
the Zhang & Kojima (2006) model prediction is consistent with 
observation at the l-cr level. 

The accretion-induced field-decay model is based on the ide- 
alized idea of dilution of polar magnetic flux due to accretion. 
All the possible instabilities are ignored. Besides, a constant 
accretion rate is assumed during the whole accretion process. 
There are, however, some speciaUties for each system, such as 
the influence of thermal and viscous instabilities of the accretion 
disk and the propeller eff'ect on the mass transfer process. The 
orbital angular momentum loss of the system and its causes may 
also have some effect on the final B and P state. The numerous 
plasma instabiUties, such as the Rayleigh- Taylor instability and 
the Kelvin-Hehnholtz instability (Ghosh & Lamb 1979a), may 
result in the penetration of the magnetosphere, prying the field 
lines aside and azimuthally wrapping the field fines by the disk 



matter (e.g. Romanova 2008; Kulkarni & Romanova et al. 2008), 
which in turn may modify the field strength evolution and then 
perturb the spin evolution. During the accretion, the accretion 
rate may change due to changes of the system and some insta- 
bilities, which can lead to the difference between the theoretical 
results and the actual values. 

It is generally believed that there are two possible ways to 
form MSPs, i.e. the standard accretion-induced field-decay and 
spin-up model and accretion-induced collapse of white dwarfs 
(AIC). Due to the conservation of magnetic flux during the col- 
lapse of white dwarfs, the MSPs formed via AIC are expected to 
have high magnetic fields and short spin periods. Although the 
number of MSPs formed via AIC is no more than 20% (Zhang 
et al. 201 1), this may also contribute to make the observed B and 
P distributions of MSPs somewhat different from that of MSPs 
formed via only the standard mechanism. 

Some selection effects should also be noted when comparing 
the model prediction with observation. For example, the spin- 
down energy loss rate, E cc ^,is related to the radio power of a 
pulsar. MSPs with relatively stronger magnetic fields and shorter 
spin periods are easier to be observed. 

Owing to the ohmic dissipation, the buried field may 
re-emerge after accretion (e.g. Young & Chanmugam 1995; 
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Fig. 4. The B and P distributions of observed MSPs (data taken 
from the ATNF pulsar catalogue). 



Bhattacharya 2008). Besides, If the magnetic field of recycled 
pulsars arrive at the bottom value at this stage and the accretion 
has not yet ended, the accumulation of material may spin up the 
recycled pulsar further. However, the magnetic field of neutron 
stars will not decay. The magnetic flux carried by the plasma 
accreted onto the neutron star surface may increase the neutron 
star surface field strength. 

We also plot the B and P evolutionary scenarios during the 
accretion process in Fig. |6] All of these plots show that the 
bottom magnetic field strength is different for different accre- 
tion rates and that the minimum period is insensitive to the ac- 
cretion rate at the end of the accretion phase. After accreted 1 
Mq, the spin period reaches about 1 ms. If the there is enough 
mass that can be accreted (e.g. 1.2 Mq), the spin period can be 
shorter than one millisecond, forming a submillisecond pulsar. 
It is claimed that the maximum accreted mass is 0.8 Mq un- 
der the assumption of 1 Mq companion star van den Heuvel & 
Bitzaraki 1995a; Wijers 1997. From the recent statistics of neu- 
tron stars/LMXB (Liu, van Paradijs & van den Heuvel 2007), 
most neutron stars/LMXBs have the companion mass of 0.7 
Mq. According to the accretion-induced field-decay and spin- 
up model, the spin period can reach the value of 1.1-1.2 ms af- 
ter accreting 0.7-0.8 Mq. So far, the observed shortest period 
for MSPs is 1.4 ms (Hessels et al. 2006), and that for millisec- 
ond X-ray pulsars is 1.6 ms (e.g. Patruno 2010). Those two spin 
frequencies, 716 Hz and 620 Hz, are both lower than than the 
believed break-up spin frequency ~ 1000 Hz (e.g. Lattimer & 
Prakash 2004). Models to explain the lack of submillisecond 
pulsars have been proposed, including angular momentum loss 
due to gravitational radiation and the magnetic spin equilibrium. 
A significant quadrupole moment may exist due to some oscil- 
lation modes or the so-called 'crustal mountains' and 'magnetic 
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Fig. 5. The computed B and P distributions of recycled MSPs 
with spin periods shorter than 20 ms. 

deformation', which leads to the emission of gravitational waves 
and the loss of angular momentum, in particularly because of 
its strong dependence on the spin frequency to the 5th power 
(e.g. Bildsten 1998; Haskefl & Pati'uno 201 1; Patruno, Haskell & 
D ' Angelo 2011). On the other hand, however, the magnetic brak- 
ing resulting from the action of stellar wind indirectly carries 
away the angular momentum (Rappaport, Verbunt & Joss 1983), 
which is more efficient than gravitational radiation by about 2 
orders of magnitude in some close systems (Kalogera, Kolb & 
King 1998). Furthermore, the magnetic spin equilibrium set by 
disk/magnetosphere coupling seems successful in explaining the 
lack of submillisecond pulsars (D' Angelo & Spruit 201 1 ; Kajava 
et al. 2011). In supplement to those efforts, the results of the 
accretion-induced field-decay and spin-up model as presented in 
this paper suggest that the achievable minimum spin period due 
to the recycling process may depend mainly on the amount of 
mass available for accretion, before the limiting shortest period 
set by the magnetic spin equilibrium is reached. 
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